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We study the dynamo threshold of a helical flow made of a mean (stationary) 
plus a fluctuating part. Two flow geometries are studied, either (i) solid body or 
(ii) smooth. Two well-known resonant dynamo conditions, elaborated for stationary 
helical flows in the limit of large magnetic Reynolds numbers, are tested against lower 
magnetic Reynolds numbers and for fluctuating flows (zero mean) . For a flow made 
of a mean plus a fluctuating part the dynamo threshold depends on the frequency 
and the strength of the fluctuation. The resonant dynamo conditions applied on the 
fluctuating (resp. mean) part seems to be a good diagnostic to predict the existence 
of a dynamo threshold when the fluctuation level is high (resp. low). 

PACS numbers: 47.65. +a 
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I. INTRODUCTION 

In the context of recent dynamo experiments [1-3], an important question is to identify the 
relevant physical parameters which control the dynamo threshold and eventually minimize 
it. In addition to the parameters usually considered, like the geometry of the mean flow 
[4, 5] or the magnetic boundary conditions [6, 7], the turbulent fluctuations of the flow seem 
to have an important influence on the dynamo threshold [8-11]. Some recent experimental 
results [12, 13] suggest that the large spatial scales of these fluctuations could play a decisive 
role. 

In this paper we consider a flow of large spatial scale, fluctuating periodically in time, 
such that its geometry at some given time is helical. Such helical flows have been identified 
to produce dynamo action [14, 15]. Their efficiency has been studied in the context of fast 
dynamo theory [16-21] and they have led to the realization of several dynamo experiments 
[3, 22-24]. 

The dynamo mechanism of a helical dynamo is of stretch-diffuse type. The radial com- 
ponent B r of the magnetic field is stretched to produce a helical field (0, B e ,B z ), where 
(r, 9, z) are the cylindrical coordinates. The magnetic diffusion of the azimuthal component 
Be produces some radial component B r due to the cylindrical geometry of the problem [17]. 
In this paper we shall consider two cases, depending on the type of flow shear necessary for 
the B r stretching. 

In case (i) the helical flow is solid body for r < 1 and at rest for r > 1 (the same 
conductivity is assumed in both domains). The flow shear is then infinite and localized at 
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the discontinuity surface r — 1. Gilbert [17] has shown that this dynamo is fast (positive 
growth rate in the limit of large magnetic Reynolds number) and thus very efficient to 
generate a helical magnetic field of same pitch as the flow. In case (ii) the helical flow is 
continuous, and equal to zero for r > 1. The flow shear is then finite at any point. Gilbert 
[17] has shown that such a smooth helical flow is a slow dynamo and that the dynamo action 
is localized at a resonant layer r = r such that < r < 1. Contrary to case (i), having a 
conducting external medium is here not necessary. 

In both cases some resonant conditions leading to dynamo action have been derived [16- 
18, 20, 21]. Such resonant conditions can be achieved by choosing an appropriate geometry of 
the helical flow, like changing its geometrical pitch. They have been derived for a stationary 
flow U(r, 9, z) and can be generalized to a time-dependent flow of the form U(r, 9, z) • f(t) 
where fit) is a periodic function of time. Now taking a flow composed of a mean part U plus 
a fluctuating part U • f(t), we expect the dynamo threshold to depend on the geometry of 
each part of the flow accordingly to the resonant condition of each of them and to the ratio 
of the intensities |U|/|U|. However we shall see that in some cases even a small intensity of 
the fluctuating part may have a drastic influence. The results also depend on the frequency 
of /(f). 

The Ponomarenko dynamo (case(i)) fluctuating periodically in time and with a fluctua- 
tion of infinitesimal magnitude had already been the object of a perturbative approach [25]. 
Here we consider a fluctuation of arbitrary magnitude. Comparing our results for a small 
fluctuation magnitude with those obtained with the perturbative approach, we found signif- 
icant differences. Then we realized that there was an error in the computation of the results 
published in [25] (though the perturbative development in itself is correct). In Appendix 
VI E we give a corrigendum of these results. 
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II. MODEL 

We consider a dimensionless flow defined in cylindrical coordinates (r, 9, z) by 

1, r < 1 

U = {0,rtt{r,t),V{r,t)) ■ h{r) with h(r) = I , (1) 

0, r > 1 

corresponding to a helical flow in a cylindrical cavity which is infinite in the z-direction, the 
external medium being at rest. Each component, azimuthal and vertical, of the dimensionless 
velocity is defined as the sum of a stationary part and of a fluctuating part 

Q(r, t) = (R m + R m f{t)) £(r), V(r, t) = (R m V + R m Vf{t)) C(r) (2) 

where R m and T (resp. R m and T) are the magnetic Reynolds number and a characteristic 
pitch of the stationary (resp. fluctuating) part of the flow. In what follows we consider a 
fluctuation periodic in time, in the form f(t) = cos(u ft). Depending on the radial profiles 
of the functions £ and ( we determine two cases (i) solid body and (ii) smooth flow 

(i) e = C = l (3) 

(ii) £ = l-r, C = l~r 2 . (4) 

We note here that the magnetic Reynolds numbers are defined with the maximum angular 
velocity (either mean or fluctuating part) and the radius of the moving cylinder. Thinking 
of an experiment, it would not be sufficient to minimize the magnitude of the azimuthal 
flow. In particular if T is large (considering a steady flow for simplicity), one would have 
to spend too many megawatts in forcing the z-velocity. Therefore the reader interested in 
linking our results to experiments should bear in mind that our magnetic Reynolds number 
is not totally adequate for it. A better definition of the magnetic Reynolds number could 
be for example R m = R m \J 1 + T 2 . For a stationary flow of type (i), the minimum dynamo 
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threshold R m is obtained for T = 1.3. 

Both cases (i) and (ii) differ in the conductivity of the external medium r > 1. In case 
(i) the magnetic generation being in a cylindrical layer in the neighbourhood of r = 1, a 
conducting external medium is necessary for dynamo action. For simplicity we choose the 
same conductivity as the inner fluid. In the other hand, in case (ii) the magnetic generation 
being within the fluid, a conducting external medium is not necessary for dynamo action, 
thus we choose an isolating external medium. Though the choice of the conductivity of the 
external medium is far from being insignificant for a dynamo experiment [3, 4, 6, 7], we 
expect that it does not change the overall meaning of the results given below. 
We define the magnitude ratio of the fluctuation to the mean flow by p = R m /R m . For 
p = there is no fluctuation and the dynamo threshold is given by R m . In the other hand 
for p ^> 1 the fluctuation dominates and the relevant quantity to determine the threshold 
is R m = pR m . The perturbative approach of Normand [25] correspond to /o < 1. 

The magnetic field must satisfy the induction equation 



where the dimensionless time t is given in units of the magnetic diffusion time, implying 
that the flow frequency Uf is also a dimensionless quantity. As the velocity does not depend 
on 9 nor z, each magnetic mode in 9 and z is independent from the others. Therefore we 
can look for a solution in the form 



where m and k are the azimuthal and vertical wave numbers of the field. The solenoidality 



dB 

~dt 



V x (U x B) + V 2 B. 



(5) 



B(r, t) = exp i{m9 + kz)b(r, t) 



(6) 



of the field V • B = then leads to 

h 77i 

-^ + b' r + i-b e + ikb z = 0. (7) 
r r 

With the new variables b ± = b r ± ibg, the induction equation can be written in the form 
f)h± j 

— h [k 2 + i(mQ + A;\/)/i(r)]6 ± = ±-rQ'h(r)(b + + b~) + £ ± 6 ± , (8) 

with 

C ±- 92 | 19 ^ ±1 ) 2 , (9 ) 

except in case (ii) where in the external domain r > 1, as it is non conducting, the induction 
equation takes the form 

(£± - k 2 ) b ± = 0. (10) 

At the interface r = 1, both B and the ^-component of the electric field E = VxB-UxB 
are continuous. The continuity of B r and Bg imply that of b ± . The continuity of B and (7) 
imply the continuity of b' r which, combined with the continuity of E z implies 

W;±^ + +n-i = o (ii) 

with D = d/dr and [h]\+ = /i( r =i_) — h( r=1+ ). We note that in case (ii) as fl r=1 - = 0, (11) 
implies the continuity of Db^ 1 at r = 1. 

In summary, we calculate for both cases (i) and (ii) the growth rate 

7 — ^(m, k, F, P, R m , R m , Uf) (12) 

of the kinematic dynamo problem and look for the dynamo threshold (either R m or R m ) 
such that the real part of 7 is zero. In our numerical simulations we shall take m = 1 
for it leads to the lowest dynamo threshold. 
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A. Case (i): Solid body flow 



In case (i) we set 



mQ + kV = R m fi + R m fif(t), with \i = m + kY and /x = m + kY, 



(13) 



and (8) changes into 



db ± 



+ [k 2 + i(R m JI + i? m /i/(t))/i(r)]6 ± = £ ± 6 ± . 



(14) 



dt 



For mathematical convenience, we take jl — . Then the non stationary part of the velocity 
does not occur in (14) any more. It occurs only in the expression of the boundary conditions 
(11) that can be written in the form 



Taking Ji — corresponds to a pitch of the magnetic field equal to the pitch of the 
fluctuating part of the flow —m/k = Y. In the other hand it is not necessarily equal to 
the pitch of the mean flow (except if Y — Y). In addition we shall consider two situations 
depending on whether the mean flow is zero (R m = 0) or not. The method used to solve 
the equations (14) and (15) is given in Appendix VIA. 

At this stage we can make two remarks. First, according to boundary layer theory results 
[16, 17] and for a stationary flow, in the limit of large R m the magnetic field which has the 
highest growth rate satisfies Ji w 0. This resonant condition means that the pitch of the 
magnetic field is roughly equal to the pitch of the flow. We shall see in section III A that 
this stays true even at the dynamo threshold. Though the case of a fluctuating flow of type 
U- f(t) may be more complex with possibly skin effect, the resonant condition is presumably 
analogous, writing Ji 0. This means that setting Ji — implies that if the fluctuations 



[Db ± ] 1 1 -± l -(R m + R m f(t))(b + + b-) 



r=l — 



0. 



(15) 
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are sufficiently large (p 1), dynamo action is always possible. This is indeed what will 
be found in our results. In other words, setting // = 0, we cannot tackle the situation of a 
stationary dynamo flow to which a fluctuation acting against the dynamo would be added. 
This aspect will be studied with the smooth flow (ii). 

Our second remark is about the effect of a phase lag between the azimuthal and vertical 
components of the flow fluctuation. Though we did not study the effect of an arbitrary 
phase lag we can predict the effect of an out-of-phase lag. This would correspond to take a 
negative value of T. Solving numerically the equations (14) and (15) for the stationary flow 
and m = 1, we find that dynamo action is possible only if AT < 0. For the fluctuating flow 
with zero mean, m = 1 and Jl = necessarily implies that kY = —1. Let us now consider a 
flow containing both a stationary and a fluctuating part. Setting T < necessarily implies 
that k > 0. Then for T > 0, the stationary flow is not a dynamo. Therefore in that case we 
expect the dynamo threshold to decrease for increasing p. For T < 0, together with T < 
and k > 0, it is equivalent to take T > and r > for k < and it is then covered by our 
subsequent results. 

B. Case (ii): Smooth flow 

For the case (ii) we can directly apply the resonant condition made up for a stationary 
flow [17, 18], to the case of a fluctuating flow. For given m and k, the magnetic field is 
generated in a resonant layer r = r where the magnetic field lines are aligned with the 
shear and thus minimize the magnetic field diffusion. This surface is determined by the 
following relation [17, 18] 

mn'(r ) + kV'(r ) = 0. (16) 
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The resonant condition is satisfied if the resonant surface is embedded within the fluid 

< r < 1. (17) 

As fl and V depend on time, this condition may only be satisfied at discrete times. This 
implies successive periods of growth and damping, the dynamo threshold corresponding 
to a zero mean growth rate. We can also define two distinct resonant surfaces r and r 
corresponding to the mean and fluctuating part of the flow, 

mU'(f ) + kV\f ) = 0, mQ'(r (t),t) + kV'(r (t),t) = (18) 

with appropriate definition of f2, V, f2 and V. In addition, if and V have the same time 
dependency, as in (2), then r becomes time independent. Then we can predict two different 
behaviours of the dynamo threshold versus the fluctuation rate p = R m /R m . If < r < 1 
and r > 1 then the dynamo threshold will increase with p. In this case the fluctuation is 
harmful to dynamo action. In the other hand if < Tq < 1 then the dynamo threshold will 
decrease with p. 

From the definitions (18) and for a flow defined by (1), (2) and (4) we have 

f = -(m/fc)/(2r) and r = -(m/k)/(2T). (19) 

For m = 1 and k < 0, taking T < implies r < and then the impossibility of dynamo 
action for the fluctuating part of the flow. Therefore, we expect that the addition of a 
fluctuating flow with an out-of-phase lag between its vertical and azimuthal components 
will necessarily be harmful to dynamo action. This will be confirmed numerically in section 
IIIC. 

To solve (8), (10) and (11), we used a Galerkin approximation method in which the 
trial and weighting functions are chosen in such a way that the resolution of the induction 
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equation is reduced to the conducting domain r < 1 [5]. The method of resolution is given 
in Appendix VI D. For the time resolution we used a Runge-Kutta scheme of order 4. 

III. RESULTS 
A. Stationary flow (R m = 0) 

We solve 

K 7 (m = 1, k, T, 0, R m , 0, 0) = (20) 

with k = (p — l)/r for case (i) and k = — l/(2f r) for case (ii). In figure 1 the threshold R m 
and the field frequency ^(7) are plotted versus JL (resp. f ) for case (i) (resp. (ii)), and for 
different values of T. Though we do not know how these curves asymptote, and though the 
range of ~jl (resp. f ) for which dynamo action occurs changes with T, it is likely that the 
resonant condition < 1 (resp. < r < 1) is fulfilled for the range of T corresponding 
to a dynamo experiment (r 1). In case (i) the dispersion relation (32) in Appendix VIA 
becomes JF = 0. 

B. Periodic flow with zero mean (R m = 0) 

We solve 

3?7(m = 1, k, 0, f , 0, R m , u f ) = 0. (21) 

In figure 2 the threshold R m is plotted versus Uf for both cases (i) and (ii). In both cases we 
take Jl = corresponding to k — — 1/F. For the case (ii) it implies from (19) that r = 1/2, 
meaning that the resonant surface is embedded in the fluid and then favourable to dynamo 
action. In each case (i) T = 1; 1.78 and (ii) T = 1;2, we observe two regimes, one at low 
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FIG. 1: The dynamo threshold R m (left column) and $$(7) (right column) versus (i) fj,, (ii) ro, for 
the stationary case, m = 1 and T = 0.5; 0.8; 1; 1.3; 2; 4; 10 . 

frequencies for which the threshold does not depend onwj and the other at high frequencies 
for which the threshold behaves like R m oc u 3 / 4 - 

To understand the existence of these two regimes, we pay attention to the time evolution of 
the magnetic field for different frequencies Uf. In figure 3, the time evolution of b~ (real and 
imaginary parts) for case (ii) V — 1 (case (c) in figure 2) is plotted for several frequencies 
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FIG. 2: Dynamo threshold R m versus ujf for case (i) with Ji = and (a) V = 1.78, (b) T = 1; for 
case (ii) with f = 0.5 and (c) f = 2, (d) f = 1. 

1. Low frequency regime 

For low frequencies (uf = 1), we observe two time scales : periodic phases of growth 
and decrease of the field, with a time scale equal to the period of the flow as expected by 
Floquet's theory. In addition the field has an eigen-frequency much higher than Uf. In fact 
the slow phases of growth and decrease seem to occur every half period of the flow. This 
can be understood from the following remarks. 

First of all the growth (or decrease) of the field does not depend on the sign of the flow. 
Indeed, from (8), we show that if 6 ± (m, k) is a solution for (Q, V), then its complex conjugate 
b (m, k) is a solution for (— fi, — V). Therefore we have 6 ± (t + T/2) = 6 ± *(t) where 
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T = 2rc/ujf is the period of the flow. Now from Floquet's theory (see Appendix VIA), 
we may write b(r,t) in the form b(r, r) exp^t), with b(r, r) 27r-periodic in r = c<j/t. This 
implies that changing (Q, V) in (— Q, — V) changes the sign of ^(7). This is consistent with 
the fact that for given m and k, the direction of propagation of B changes with the direction 
of the flow. Therefore changing the sign of the flow changes the sign of propagation of the 
field but does not change the magnetic energy, neither the dynamo threshold R m which are 
then identical from one half-period of the flow to another. This means that the dynamo 
threshold does not change if we consider f(t) = |cos(u;/i)| instead of cos(u ft). It is then 
sufficient to concentrate on one half-period of the flow like for example [^y, ^y] (modulo 

TT). 

The second remark uses the fact that the flow geometry that we consider does not change in 
time (only the flow magnitude changes). For such a geometry we can calculate the dynamo 
threshold R m corresponding to the stationary case. Then coming back to the fluctuating 
flow, we understand that R m \f(t)\ > R m (resp. R m \f(t)\ < R m ) corresponds to a growing 
(resp. decreasing) phase of the field. Assuming that the dynamo threshold R m is given by 
the time average < . > of the flow magnitude leads to the following estimation for R m 

~ 7T 

Rm ~ -^Rrn (22) 

as < |cos(u;/i)| >= 2/ir. For the three cases (a), (b) and (c) in figure 2 we give in table 
III B 1 the ratio 2R m /irR m which is found to be always close to unity. In this interpretation 
of the results the frequency ujf does not appear, provided that it is sufficiently weak in order 
that the successive phases of growth and decrease have sufficient time to occur. This can 
explain why for low frequencies in figure 2 the dynamo threshold R m does not depend on 

Uf. 

Finally the frequencies u of the stationary case for T = T are also reported in table III B 1. 
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TABLE I: see in the text. 

For a geometry identical to case (c), we find, in the stationary case, UJ = 33 which indeed 
corresponds to the eigen- frequency of the field occurring in figure 3 for ojj = 1. The previous 
remarks assume that the flow frequency is sufficiently small compared to the eigen-frequency 
of the field, in order to have successive phases of growth and decrease of the field. We can 
check that the values of u given in table III B 1 are indeed reasonable estimations of the 
transition frequencies between the low and high frequency regimes in figure 2. 

2. High frequency regime 

In case (ii) and for high frequencies (figure 3, cuf = 100), the signal is made of harmonics 
without growing nor decreasing phases. We note that the eigen-frequency of the real and 
imaginary parts of b~ are different, the one being twice the other. 

In case (i), relying on the resolution of equations (14) and (15) given in Appendix VI A, we 
can show that R m oc uj 3 ^ 4 \ We also find that some double frequency as found in figure 3 for 
the case (ii) can emerge from an approximate 3x3 matrix system. As these developments 
necessitate the notations introduced in Appendix VIA, they are postponed in Appendix 
VI B. 
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3. Further comments about the ability for fluctuating flows to sustain dynamo action 

We found and explained how a fluctuating flow (zero mean) can act as a dynamo. We 
also understood why the dynamo threshold for a fluctuating flow is higher than that for 
a stationary flow with the same geometry. It is because the time-average of the velocity 
norm of the fluctuating flow is on the mean lower than that of the stationary flow. This can 
be compensated with other definitions of the magnetic Reynolds number. Our definition is 
based on max t \ fl(r, t) \ . An other definition based on < \ ft(r, t) \ > t would exactly compensate 
the difference. 

Recently a controversy appeared about the difficulty for a fluctuating flow (zero mean) to 
sustain dynamo action at low P m [26] , whereas a mean flow (non-zero time average) exhibits 
a finite threshold at low P m [10, 27] (the magnetic Prandtl number, P m , being defined as 
the ratio of the viscosity to the diffusivity of the fluid). This issue is important not only 
for dynamo experiments but also for natural objects like the Earth's inner-core or the solar 
convective zone in which the electro-conducting fluid is characterized by a low P m . Though 
we did not study this problem, our results suggest that the dynamo threshold should not 
be much different between fluctuating and mean flows, provided an appropriate definition 
of the magnetic Reynolds number is taken. In that case why does it seem so difficult to 
sustain dynamo action at low P rn for a fluctuating flow [26] whereas it seems much easier 
for a mean flow [10, 27] ? 

In the simulations with a mean flow [10, 27], two dynamo regimes have been found, 
the one with a threshold much lower than the other. In the lowest threshold regime, the 
magnetic field is generated at some infinite scale in two directions [27]. There is then an 
infinite scale separation between the magnetic and the velocity field, and the dynamo action 
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is probably of mean-field type and might be understood in terms of a-effect, /3-effect, etc. 
In that case, removing the periodic boundary conditions would cancel the scale separation 
and imply the loss of the dynamo action. In the highest threshold regime, the magnetic 
field is generated at a scale similar to the flow scale, the periodic boundary conditions are 
forgotten and the dynamo action can not be understood in terms of an ct-effect anymore. 
In order to compare the mean flow results [27] with those for a fluctuating flow [26] we have 
to consider only the highest threshold regime in [27], the lowest one relying on mean-field 
dynamo processes due to the periodic boundary conditions and which are absent in the 
fluctuating flow calculations [26]. 

Now when comparing the threshold of the highest threshold regime for a mean flow with 
the threshold obtained for a fluctuating flow and with appropriate definitions of R m , a strong 
difference remains at low P m [28]. A speculation made by Schekochihin et al. [29] is that 
the highest threshold regime obtained for the mean flow at low P m [27] would correspond in 
fact to the large P m results for the fluctuating flow in [28]. Their arguments rely on the fact 
that the mean flow in [27] is peaked at large scale and so spatially smooth for the generated 
magnetic field. It would then belongs to the same class as the large- P m fluctuation dynamo. 
Both dynamo thresholds are found to be similar indeed, and thus the discrepancy vanishes. 

Though the helical flow that we consider here is noticeably different (no chaotic trajec- 
tories) it may have some consistency with the simulations at large P m mentioned above and 
at least supports the speculation by Schekochihin et al. [29]. 

C. Periodic flow with non zero mean 

We are now interested in the case where both R m ^ and R m ^ 0. The flow is then the 
sum of a non zero mean part and a fluctuating part. We have considered two approaches 
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depending on which part of the flow geometry is fixed, either the mean or the fluctuating 
part. 

i. r = 1 

Here we fix F — 1, m — 1 and k — — 1 and vary T, p and cj/ for the case (ii). Then we 
solve the equation 

gfc 7 (m = l, fc = -l, r = l, f = l/2?o, i? m , R m = pR m , u f ) = Q (23) 

to plot i? m as a function of p in figure 4 for values of r and uif. From (19) we have f = 1/2 
which corresponds to a mean flow geometry with a dynamo threshold about 100. The curves 
are plotted for several values of T leading to values of r not necessarily between and 1. 
We consider two fluctuation frequencies ujf — 1 and uof = 50. We find that the dynamo 
threshold R m increases asymptotically with p unless the resonant condition < r < 1 is 
satisfied, here the curves (a), (b) and (e). For these three curves we checked that in the 
limit of large p, R m = 0(p~ 1 ). For r = 1/4 (curve (e)) and for p as 1 we do not know if a 
dynamo threshold exists. 



2. r = i 

Here we fix F — 1, m — 1 and k = — 1 and vary F, p and a;/. We then solve the equation 

sft 7 ( m = i ? k = -l, T, f, i? m , R m = pR m , uj f ) = (24) 

with r = 1 — /7 in case (i) and T = l/2r in case (ii). In figure 5, R m is plotted versus p 
for values of pZ (resp. f ) in case (i) (resp. (ii)) and Uf. Taking r = 1, m = 1 and k = — 1 
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FIG. 4: Dynamo threshold R m versus p for case (ii), for two frequencies LOf = 50 and loj = 1 and 
ro = 1/2 (r = 1, m = 1, k = —1). The different curves correspond to fo = (a) 1/2; (b) 2/3; (c) 1; 



(d) oo; (e) 1/4; (f) -1; (g) -1/2 (r = (a) 1; (b) 0.75; (c) 0.5; (d) 0; (e) 2; (f) -0.5; (g) -1). 



implies Ji = in case (i) and r = 0.5 in case (ii). In both cases (i) and (ii) the fluctuating 
part of the flow satisfies the resonant condition for which dynamo action is possible. This 
implies that R m should scale as 0(p^ 1 ) provided that p is sufficiently large. In each case 
we consider two flow frequencies, uif — 1 and Uf = 10 for case (i), cuf — 1 and Uf = 50 
for case (ii). The curves are plotted for different values of T corresponding to < 1 for 
case (i) and < r < 1 for case (ii). For large p we checked that R m = 0(p^ 1 ). The main 
differences between the curves is that R m versus p may decrease monotonically or not. In 
particular in case (i) for Jl = 0.4, R m decreases by 40% when p goes from to 1 showing 
that even a small fluctuation can strongly decrease the dynamo threshold. In most of the 
curves there is a bump for p around unity showing a strong increase of the threshold before 
the final decrease at larger p. 
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FIG. 5: The dynamo threshold R m versus p for k = —1, m = 1 and r = 1 (// = in case (i) and 
?o = 0.5 in case (ii)) and ujf = 1, 10 or 50. The labels correspond to ~p in case (i) and vq in case 
(ii). 

3. r = f = i 

Here we fix T = T = 1, m = 1 and = — 1 for the case (ii) and vary ujf and p. We then 
solve the equation 



K 7 (m=l, fc = -l, r = l, T = l, i? m , R m = pR m , u f ) = 



(25) 
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FIG. 6: Dynamo threshold R m versus p for tq = tq = 0.5 (T = T = 1). The labels correspond to 
different values of ujf. The eigen-frequency for p = is ZJ = 33. 

to plot R m versus p in figure 6 for various frequencies Uf. Taking r = r = l,m = l and 
k — — 1 implies f = r = 0.5 For p larger than 1, i? m decreases as 0(p _1 ) as mentioned 
earlier. For p smaller than unity, R m decreases versus p monotonically only if Uf is large 
enough. In fact the transition value of Uf above which R m decreases monotonically is exactly 
the field frequency a; (here a; = 33) corresponding to p = . This shows that a fluctuation 
of small intensity (p < 1) helps the dynamo action only if its frequency is sufficiently high. 
This is shown in Appendix VI C for the case (i). Though, the frequency above which a small 
fluctuation intensity helps the dynamo may be much larger than uj. For example in case 
(i) for Jl = 0.4 and Jl = represented in figure 5, we have Zu = 0.51. For Uf — 1 small 
fluctuation helps, for Uf = 10 they do not help, and for higher frequencies they help again. 
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IV. DISCUSSION 

In this paper we studied the modification of the dynamo threshold of a stationary helical 
flow by the addition of a large scale helical fluctuation. We extended a previous asymptotic 
study [25] to the case of a fluctuation of arbitrary intensity (controlled by the parameter 
p). We knew from previous studies [17, 18] that the dynamo efficiency of a helical flow is 
characterized by some resonant condition at large R m . First we verified numerically that 
such resonant condition holds at lower R m corresponding to the dynamo threshold, for both 
a stationary and a fluctuating (no mean) helical flow. Then for a helical flow made of a mean 
part plus a fluctuating part we showed that, in the asymptotic cases p < 1 (dominating 
mean) and p 1 (dominating fluctuation), it is naturally the resonant condition of the 
mean (for the first case) or the fluctuating (for the second case) part of the flow which 
governs the dynamo efficiency and then the dynamo threshold. In between, for p of order 
unity and if the resonant condition of each flow part (mean and fluctuating) is satisfied, the 
threshold first increases with p before reaching an asymptotic behaviour in 0(p _1 ). However 
there is no systematic behaviour as depicted in figure 5 (i) for cuf — 1 and pZ = 0.4, in which 
a threshold decrease of 40% is obtained between p = and p — 1. If the fluctuation part 
of the flow does not satisfy the resonant condition, then the dynamo threshold increases 
drastically with p. 

Contrary to the case of a cellular flow [11], there is no systematic effect of the phase lag 
between the different components of the helical flow. For the helical flow geometry it may 
imply an increase or a decrease of the dynamo threshold, depending how it changes the 
resonant condition mentioned above. 

There is some similarity between our results and those obtained for a noisy (instead of 
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periodic) fluctuation [8]. In particular in [8] it was found that increasing the noise level 
the threshold first increases due to geometrical effects of the magnetic field lines and then 
decreases at larger noise. This could explain why at p 1 we generally obtain a maximum 
of the dynamo threshold. 

Finally these results show that the optimization of a dynamo experiment depends not only 
on the mean part of the flow but also on its non stationary large scale part. If the fluctuation 
is not optimized then the threshold may increase drastically with (even small) p, ruling out 
any hope of producing dynamo action. In addition, even if the fluctuation is optimized 
(resonant condition satisfied by the fluctuation), our results suggest that there is generally 
some increase of the dynamo threshold with p when p < 1. If the geometry of the fluctuation 
is identical to that of the mean part of the flow, there can be some slight decrease of the 
threshold at high frequencies but this decrease is rather small. When p > 1 the dynamo 
threshold decreases as 0(p~ l ) which at first sight seems interesting. However we have to 
keep in mind that as soon as p > 1 the driving power spent to maintain the fluctuation 
is larger than that to maintain the mean flow. Then the relevant dynamo threshold is not 
R m any more, but R m = pR m instead. In addition monitoring large scale fluctuations in an 
experiment may not always be possible, especially if they occur from flow destabilisation. 
In that case it is better to try cancelling them as was done in the VKS experiment in which 
an azimuthal belt has been added [30]. 
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VI. APPENDIX 

A. Resolution of equations (14) and (15) for the case (i): solid body flow 

As f(t) is time-periodic of period 2ir/ujf, we look for b(r, t) in the form b(r, r) exp(7t) 
with b(r, r) being 27r-periodic in r = Uft. Thus we look for the functions b ± (r,r) in the 
form 

b ± (r,r) = b^(r) exp(mr) (26) 
where, from (14) and for /I = 0, the Fourier coefficients b^(r) must satisfy 

[ 7 + k 2 + i(R m Jlh{r) + nu) f )]b± = CH*. (27) 

In addition, the boundary condition (15) with f(t) = cos(t) implies 

[Db±]\- ± l -R rn (b+ + b-) r=1 ± l -R m (bti + K-i + K + i + ^+i)r=i = 0. (28) 

The solutions of (27) which are continuous at r = 1 can be written in the form 

= Cfitf, with tf=l , (29) 

[K^M/K^sn), r>l 

with 

q 2 = k 2 + 7 + i(Rmp + nuf), s 2 n = k 2 + 7 + inu f . (30) 
Substituting (29) in (28), we obtain the following system 

CX ± i^(C: + C~) ± i^{CU + + C: +1 + C" +1 )) = (31) 

with K± = qJn l l t ~ s nK'nl K t and where J n = I m±i{q n ) and = K m±1 (s n ) are 
modified Bessel functions of first and second kind. 
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The system (31) implies the following matrix dispersion relation 

TnC n ~ i^(K ~ ^n)(Cn-l + C n+1 ) = (32) 

with Cj = C+ + CJ and 

T n = 11+11- - i(R m /2)(1l+ - 1Z-). (33) 

Solving the system (32) is equivalent to setting to zero the determinant of the matrix A 
defined by 

R 

A nn — J~n 4n-i = A n n+ \ = —i—^-(1Z+ — 1Z n ) (34) 
and with all other coefficients being set to zero. 

B. High frequency regime for the periodic flow (i) with zero mean 

Following the notation of section VIA, and considering a periodic flow (i) with zero mean, 
we have JI = 0. From (30) this implies that q n = s n . Using the identity 

l^K - = f (35) 

we obtain 1Z+ = {I+K+y 1 . As R m = 0, the equation (33) becomes T n = l^n^n- Then we 
can rewrite the system (32) in the form 

C n + i^(I+K+ - I-K-){C n ^ + C n+1 ) = 0. (36) 

From the asymptotic behaviour of the Bessel functions for high arguments, we have 

a n = I-K- - I+K+ « 1/sl (37) 

For the high values of n these terms are negligible and in first approximation we keep in the 
system (36) only the terms corresponding to n — 0, ±1. This leads to a 3 x 3 matrix system 
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whose determinant is 



1 + 



(Rr, 



16 



-Qfo^-l + a l) — 0. 



(38) 



At high forcing frequencies Uf we have s±i ~ ^/w/. Together with (38), it implies 



p _ , 3/4 



(39) 



In addition, from the approximate 3x3 matrix system, the double- frequency 2u/ emerges 
for n — ±1. 

C. High frequency regime and small modulation amplitude for the periodic flow (i) 

with non zero mean 



For small amplitude modulation p « 1 the system (31) is truncated so as to keep the 
first Fourier modes n = and n — ±1. The dispersion relation : 



W(-f) «-k -)(-^-i + ^-^)=o 



(40) 



is then solved perturbatively setting R m = R + 5R and U = oj + 5u and expanding 
J r o(R m ,u)) to first order in 5R and <5u; given that jF (i? , cj ) = and with the constants 
= ±1Zq . The dispersion relation (40) becomes 



dR r 



did 



R ( 



+i 



(41) 



with [3 n = — 1Z n . The threshold and frequency shifts which behave like p 2 are written 
SR = p 2 i?2 and 5uo = p 2 oo 2 - In the left-hand-side of (41) the partial derivatives are given by 



duo 
dR m 



1 

"Co 
1 

Co 



(cty 



duj 

ORrr, 



(Co) 



duj 
> dKo 
dRm 



-Co 



(42) 
(43) 
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One can show that the partial derivatives of TZq are related to integrals calculated in [25] 
through the relations 

° " i I (^?rdr ^=^Jl (*J) 2 rdr (44) 

m JO 



du L y v/ OR 



In the following we shall focus on the case p = and we introduce the notations 

i(\M. I ih). ^ = -iCotgi + ig2), ^ + ^ = X + iY (45) 



dR m uyJ± J * n du '~^ ax T-x T. 



+i 



Solutions of (41) are 

„ fRv\ 2 Xg 1+ Yg 2 ( RoV Xfi+Yfr 

R " = [t) " 2 = [t) (46) 

recovering results similar to those obtained in [25] using a different approach. 

We have in mind that for some values of ujf resonance can occur. An oscillating system 
forced at a resonant frequency is prone to instability and a large negative threshold shift 
is expected. However, inspection of (46) reveals no clear relation between the sign of R2 
and the forcing frequency which appears in the quantities X and Y . We only know that 
/i<72 — fi9i < 0, since near the critical point (5R = R m — R ) the denominator in (46) is 
proportional to the growth rate of the dynamo driven by a steady flow. When p = 0, we 
shall consider Eq. (41) for an imposed 5R and complex values of 5uj = oj\ + io\ where lo\ is 
the frequency shift and Sft(7) = —<J\ is the growth rate, given by 

x D /i<72 — hgi /._s 
ai = 5R n 2 + n 2 ( 47 ) 

i/l "T i/2 

Above the dynamo threshold, (5R > 0) the field is amplified (^(7) > 0) thus U\ < and 
fi92 - /201 < 0. 

In the high frequency limit (ujf » Uo) expressions for X and Y can be derived explicitly 
using the asymptotic behaviour of the Bessel functions for large arguments. For pZ — with 



28 

q±i = s±i (c<j/±c<Jo) 1 ^ 2 (l±0/v / 2 and using the asymptotic behaviour : P±i/J 7 ±i — > (s±i)~ 3 
one gets 

X + ^ = -^7 3/2 (l-^^) (48) 

2cjj 

When /7 = 0, we have also fx = 1/2 and f 2 = 0, leading to the expression for R 2 : 

it> 2 « -^T 3 ' 2 - ^ (49) 

For the wave numbers m = —k = 1, numerical calculations of g\ and g 2 which only depend 
on the critical parameters R and u give g±/g 2 = 1.626, and thus R 2 < when — > 00. 

When 77 = 0, there are several reasons to consider the particular value of the forcing : ujf — 
2u)q. One of them is that for Hill or Mathieu equations it is a resonant frequency. Moreover, in 
the present problem it leads to simplified calculations. In particular the asymptotic behavior 
of p n / ' T n can still be used for n = +1 since ujf + u is large, while the approximation is no 
longer valid for n — — 1. Nevertheless, the mode n = — 1 is remarkable since it corresponds 
to s_i = Sq from which it follows that = (3q and = —iR f3^. Finally one gets the 
exact result : P-\jT-\ = i/Ro, which leads to 

3/2 



X + iY^^- + \ with s\^-2(^f] {l-i 
Ro s 1 \ 2 J 



) (50) 



For the values R = 20.82 and c^o = 4.35 corresponding to Fig. 8 (f) one gets X = — 1.5 x 10 2 
and Y = 3.3 x 10~ 2 . The threshold shift is 

d2 

R 2 w -^(1.62X + y) = 0.48 (51) 
8 

showing that the sign of i? 2 changes when a;/ decreases from infinity to 2u . This result 
is in qualitative agreement with the exact results reported in Fig. 8 (f) where R 2 = for 
Uf = 8.3 ~ 2u . When the forcing frequency is exactly twice the eigen-frequency u we 
had rather expected a large negative value of R 2 on the basis it is a resonant condition 
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for ordinary differential system under temporal modulation. In Fig. 8 (f) the maximum 
negative value of R 2 occurs for uj ~ 4cj which cannot be explained by simple arguments. 

For fx ^ 0, we have not been able to find resonant conditions like : nujf + mw = 
(n, m integers) between u>f and u such that Uf would be associated to a special behavior 
of the threshold shift. Contrary to the Hill equation, the induction equation is a partial 
differential equation with the consequence that the spatial and temporal properties of the 
dynamo are not independant. The wave numbers k and m are linked to the frequencies u 
and ujf through q± n and s± n which appears as arguments of Bessel functions having rules of 
composition less trivial than trigonometric functions. Exhibiting resonant conditions implies 
to find relationship between q± n , s± n and q , s for specific values of ujf. We have shown 
above for JI — that a relation of complex conjugaison exists for n = —1 when ujf — 2u 
but we have not yet found how to generalize to other values of JI and have left this part for 
a future work. 



D. Resolution of case (ii): smooth flow 

We define the trial functions vpf = K m ±i(k)J m ±i(ajr) / J m ±i{®-j) where the atj are the 
roots of the equation 

\ K m+ i(k) _ Km-^k) 

1 \_Jm+l( a j) Jm-l( a j) 

and where J and K are respectively the Bessel functions of first kind and the modified Bessel 
functions of second kind. For r < 1, we look for solutions in the form 

N 

p = Y,»wt( a > r ) ( 53 ) 



+ 2k 



K m {k) 



o, 



(52) 
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where N defines the degree of truncature. For r > 1 the solutions of (10) are thus of the 
form 

N 

b ± = K m±1 {kr)Y,V{t) (54) 

3=1 

and, from (52) these solutions satisfy the conditions (11) at the interface r = 1. To determine 
the functions bj(t) it is sufficient to solve the induction equation (8) for r < 1. For that 
we replace the expression of b^ 1 given by (53) into the induction equation (8), in order to 
determine the residual 

N . N 

R± = J2 + ^ + a ) + l[m£l + w )\ h i) T \rQ! h 3^t + ^j)- ( 55 ) 

3=1 ~ 3=1 

Then we solve the following system 

/ R + <ptrdr + / RTfcrdr = 0, i = 1, • • • , N (56) 
Jo Jo 

where the weighting functions are defined by (fif = J m ±i( a j r )/ Jm±i( a j)- Using the orthog- 
onality relation 

/ <pjrfrdr+ [ <P~i)-rdr = dijdj (57) 
Jo Jo 

with 

Gu = K j™ +1 j\ \ [ Jm+Mir)rdr + [ J 2 m _ x {a,<r)rdr, (58) 

J m+l\ a i) JO J m-l\ a i) JO 

we write the system (56) in the following matrix form 

X = MX with X = (6i, ■■■ ,b N ) (59) 

with 

Mlj = Si^k 2 + a)) + / (0+^ + + <j>;i/>j)(mSl + W)rdr (60) 
- ^ / V - W" + V-D^dr. (61) 

^^ii JO 

The numerical resolution of this system is done with a fourth order Runge-Kutta time-step 
scheme. We took a white noise as an initial condition for the bj. 
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E. Corrigendum of Normand (2003) results [25] 



For very small values of the fluctuation rate p and for an infinite shear (case (i)), a 
comparison can be made between the results obtained for Jl — by the method based on 
Floquet theory (see Appendix VIA) and those obtained by a perturbative approach [25] 
which consists in expanding R m and the frequency S(p) in powers of p according to 



where Rq and uoq are respectively the critical values of the Reynolds number and the fre- 
quency in the case of a stationary flow. At the leading order it appears that : lo\ — R\ — 0. 
At the next order, the expressions of R 2 and uj 2 are given in [25], however their numerical 
values are not correct due to an error in their computation. 

After correction, the new values of R 2 and oj 2 are given in Figure 7 for the set of parameters 
considered in [25]: m — 1, k — —0.56 and T = 1 (p = 0.44). The different curves correspond 
to values of T which are not necessarily the same than those taken in [25]. For \p\ sufficiently 
small (curves e and f), R 2 changes its sign twice versus the forcing frequency. We find that 
R 2 is negative for low and high frequencies implying a dynamo threshold smaller than the 
one for the stationary flow. At intermediate frequencies R 2 is positive with a maximum 
value, implying a dynamo threshold larger than the one for the stationary flow. For larger 
values of \p\ (curves a, b, c and d) R 2 is positive for all forcing frequencies with a maximum 
value at a low frequency which increases with \]x\. This implies that for sufficiently large 
the dynamo threshold is larger than the one obtained for the stationary flow as was already 
mentioned in section III A. For Y = 1.78 we have p, — 0. In this case, we have checked that 
the values of R 2 and u 2 are in good agreement with the values of R m and uj obtained by the 



$t(p) = UJ + pUJi + p 2 UJ 2 + .... 



(62) 



R m = R + pR t + p 2 R 2 + .... 



(63) 
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FIG. 7: Results obtained by the perturbative approach for m = 1, k = —0.56 and ll = 0.44 (r = 1). 
The dynamo threshold i?2 is plotted versus ujj for several values of ju = (a) 1.56; (b) 1.28; (c) 1; 
(d) 0.72; (e) 0.44; (f) 0.16 and (g) 0. 




FIG. 8: Same as Figure 7 but for m = 1, k = — 1 and p = (r = 1). The labels correspond to Jl = 
(a) 1.5; (b) 1.25; (c) 1; (d) 0.5; (e) -0.5 and (f) 0. 

method of Appendix VIA, provided p < 0.1. 

For completeness we have also calculated the values of R 2 and o; 2 for m — 1, k — — 1 and 
r = 1, (/7 = 0) as considered in the body of the paper. The results are plotted in Figure 
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8. Qualitatively the results are in good agreement with those of Figure 7. For jj, — again, 
we have checked that the values of R 2 and u 2 are in good agreement with the values of R m 
and ZJ obtained by the method of Appendix VIA, provided p < 0.1. For higher values of 
the modulation amplitude p the relative difference between the results obtained by the two 
methods can reach 10% on R 2 for p = 0.4. 

Finally it must be noticed that our parameters T and Uf are strictly equivalent to respec- 
tively E\ and a in [25]. 
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TABLE II: see in the text. 



